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We describe analytical and numerical results on the statistical properties of complex eigenvalues 
and the corresponding non- orthogonal eigenvectors for non-Hermitian random matrices modeling 
one-channel quantum-chaotic scattering in systems with broken time-reversal invariance. 
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The statistical properties of non-orthogonal eigenvec- 
tors of large non-self adjoint random matrices have re- 
cently been characterised in Refs. P^,^,Hj5[ 

Correlations of non-orthogonal eigenvectors are ex- 
pected to determine dynamical properties of classical ran- 
dom systems described by njon-selfadjoint operators, such 
as Fokker-Planck operatorsQ for example; they also play 
an important role in quantum systems: in Ref. |^ it was 
observed that the statistics of non-orthogonal eigenvec- 
tors determines the properties of random lasing media. 
This has led to an increased interest in eigenvector statis- 
tics in non-self adjoint random matrix ensembles (see also 
Ref. 0). 

In a model for quantum-chaotic scattering, the com- 
plex eigenvalues £k ,k = 1, ...,N of a random N x N 
non-Hermitian matrix (the so-called "effective Hamil- 
tonian") Hn — H — iV are used to describe generic 
statistical properties of resonances in quantum chaotic 
scattering (see Ref. || and references therein): for sys- 
tems with broken time-reversal invariance (anti-unitary 
symmetry) , the matrices H are random Pi,x N matri- 
ces from the Gaussian Unitary Ensemblecl with joint 
probability density P{H)dH cx exp[-(7V/2)Tri72] rfij. 
In the limit of large N , the mean eigenvalue density 
y{E) for such matrices is given by the semicircular law 
v{E) = (27r)-V4-£;2 for \E\ < 2 (and zero other- 
wise). The corresponding mean spacing between neig- 
bouring eigenvalues around the point E in the spectrum 
is given by A{E) = l/[Nv{E)]. 

The Hermitian matrices H describe the energy-level 
statistics of the closed counterpart of the scattering sys- 
tem; the Hermitian NxN matrix F > models the cou- 
pling of the system to scattering continua via M = 1, 2, ... 
open channels. It has rank M < N. For our purposes it 
can be chosen diagonal: F — diag(7i, 72, 7a/, 0, 0). 
The constants < 7c < 00 parametrise the strength of 
the coupling to the scattering continua via a given chan- 
nel c = l,...,Af. Here 7c = corresponds to a closed 
channel c, and 7c = 1 describes the so-called perfectly 
coupled channel. Empirical situations correspond to the 
regime of large N, with M fixed and M ^ N. Then the 
widths Ffe = 2 Im £k are of the same order 1 /N as the 
mean spacing A{E) between the positions of the neigh- 
bouring resonances along the real energy axis. In this 
regime, the resonances may partly or considerably over- 



lap and first-order perturbation theory valid for small res- 
onance overlaps brea|ks,|dflwn. Similarly, self-consistent 
perturbation schemealralj^El assuming many channels and 
strongly overlapping resonances are inapplicable. 

A detailed analytical understanding of the statistical 
properties of the resonances in the regime of partial over- 
lap has recently been achieved for ijti|e-|Case of systems 
with broken time reversal invarianceolla. These results, 
based on the random matrix approach, are expected to be 
applicable to a broad class of quantum-chaotic systems. 
Indeed, the distribution of the widths F^ derived in Ref. ^ 
is in good agreement with available numerical data. .for 
quite diverse models of quantum chaotic scatteringljuj. 

Much less is known on properties of non-orthogonal 
eigenvectors. Let \Rk) and {Lk\ denote the right and the 
left eigenvectors of the matrix Ti. corresponding to the 
eigenvalue £k = Ek — iYk = Ek — iTk/2, 



'H\Rk) — £k\Rk) 
n^\Lk) = £l\Lk) 



{Lk\n ={Lk\£k 
{Rk\n^ ^ {Rk\£*k 



(1) 



where the symbols f and * stand for Hermitian con- 
jugation and complex conjugation, respectively. Ex- 
cept for a set of measure zero, the eigenvalues are non- 
degenerate. In this case the eigenvectors form a com- 
plete, bi-orthogonal set. They can be normalised to sat- 
isfy {Lk\Ri) = Ski ■ The most natural way to charac- 
terise the non-orthogonality of eigenvectors is to consider 
statistics of the overlap matrix Oki = {Lk\Li){Ri\Rk) ■ 
This matrix features in two-point correlation functions 
in non-Hermitian systems, e.g. in description of the par- 
ticle escape from the scattering region ( "norm leakage" , 
see Ref. 

Following Ref. |^, consider two correlation functions: a 
diagonal one 



0{£)^{];^Y.'^nn5{£~£k 



(2) 



and an off-diagonal one 



0{£i,£2) = (^E^-- 5{£^-£n) 5{£2-£„,) ) . (3) 



Here (• • ■)'Hn stands for an ensemble average over T-Ln- 
The correlation functions (0,0) characterise the average 
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non-orthogonality of eigenvectors corresponding to reso- 
nances whose positions in the complex plane are close to 
the complex energies £, and Here stands for 

a two-dimensional 5— function of the complex variable £ . 

In the context of lasing media, the diagonal correlator 
characterises average excess noise factors (Petermann 
factors), and the off-diagonal correlator (0) describes 
average cross correlations between themial or quantum 
noise emitted into different eigenmodeaiJ. Note that for 
any ensemble with orthogonal eigenvectors and complex 
eigenvalues E (for normal matrices), 0(5) is equal to 
the mean density of complex eigenvalues, and the off- 
diagonal correlator vanishes: 0(£i,£2) = 0. 

Both diagonal and off-diagonal eigenvector correlators 
were introduced and calculated for the case of Ginibre's 
ensemble of non-Hermitian matrices in Ref. |l|. For the en- 
semble Hat pertinent to chaotic scattering, both types of 
eigenvector correlators were found recently for the regime 
of very strongly overlapping resonances wbeii widths typ- 
ically much exceed the mean separationB'El. Physically 
this regime corresponds to a situation where the scat- 
tering system is coupled to the coatinuum via a large 
number M ^ 1 of open channelaa. In this casfii-the 
self-consistent Born approximation is adequateEl'B'ffl, a 
perturbative approximation valid in the limit of large iV, 
large M, and |£i — £2! 7^ 0, provided E\,£i are well in- 
side the support of the spectrum. A non-perturbative 
expression for the diagonal correlator 0(z) valid for any 
number of open channels was obtained in Ref. |^ by em- 
ploying a heuristic analytic continuation procedure. For 
the case of the resonance widths, this heur^tic scheme 
is known to reproduce the exact expressionQ. It is thus 
natural to expect that this procedure is adequate in the 
case of eigenvector statistics, too, although this remains 
to be proven. 

No non-perturbative results for the off-diagonal eigen- 
value correlator 0{£i,£2) have so far been reported, to 
the best of our knowledge. 

In the present paper we provide exact non-perturbative 
expressions for both diagonal and off-diagonal eigenvec- 
tor correlators valid for the case of a system with broken 
time-reversal invariance (anti-unitary symmetry) cou- 
pled to continuum via a single open channel (M = 1) 
with coupling strength 7. The single-channel case de- 
scribes pure resonant chaotic reflection. This case is 
more amenable to analytical treatment than a general 
case (M > 1), combining both reflection and transmis- 
sion phenomena. Understanding the single-channel case 
should be considered as a useful step towards a more 
complete pictureEj. 

Our result for the diagonal correlator is 
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sinh (27ry/A) 1 
27rr/A J 



(4) 



where E ^ E - iY , v = v{E) , A = 1\{E) and 
.9 = (7 + 7~^)/(27ri^) is the effective (renormalised) cou- 
pling strength. The result for 0{£) agrees with one 
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FIG. 1: Numerical (□) and analytical results (solid line) for 

(a) d{£) and 0{£) as a function of Y for A'' = 32, and 7 = 0.9; 

(b) {n{L^,Ly)) and Oi{L^,Ly) for N = 32, = 0.1, and 
7 = 0.9; (c) 02iL^,Ly) for iV = 32, = 0.1, and 7 = 0.9. 



reported in Ref. ^ confirming the validity of the ana- 
lytical continuation scheme used there. For the sake 
of comparison we present here also the expression for 
the single-channel resonance density xlefined as d{£) = 
(iV-i Efc ^ - and given byB: 



-2.,y/A sinh(2^y/A) \ 
27rr/A / ■ 



We have compared these analytical expressions, valid 
in the limit iV — > 00, with direct numerical diagonali- 
sations of finite-dimensional matrices Hn , see Fig. ^a) . 
This is of interest since empirically, the ensemble average 
(•••)'«« is usually replaced by an energy average over sev- 
eral spectral windows, each of which may typically con- 
tain of the order of 10 or 100 resonances, corresponding 
to a finite value of N. We observe that the analytical 
results describe the numerical data well, except for small 
deviations at large values of Y. Numerically it is easier to 
compute smoothed averages, such as the mean number 
of eigenvalues {n(Lx, Ly)) inside a rectangular domain 



A 



<Re£ < 
< lm£ < L„ 



(6) 



in the complex plane. This quantity can be obtained from 
the mean density d(£) by integration over A. Similarly 
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FIG. 2: Numerical (□) and analytical results (solid line) for 
the distribution of the narrowest resonance for l3 — 2, and 
7 = 0.1, W = 0.2, and N = 128. Here x = ngnV/A. 



one can define the function Oi{Lx, Ly) as the integral of 
the diagonal correlator 0{£) over the same domain, ob- 
taining Oi{Lx,Ly) = {N~^ Se^GA Okk)- Numerical ver- 
sus analytical results for these two quantities are plotted 
in Fig. |l|(b). 

For the off-diagonal correlator 0{£i,£2) we obtain 
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7ry2(Ai-A2)/A 



(7) 



,7rYi(Ai-A2)/A 



-7rYi(Ai-A2)/A 



where Re £1^2 = ^'1,2 = E ^ Q, and it is assumed that 
57 ~ A. We have also calculated the corresponding 
smoothed average 02{Lx,Ly) = {N~'^J2£^^£„eA^"in), 
by integrating £1 and £2 in over the domain A. In 
Fig. 0(c) we compare this result (valid in the limit of 
A'^ 00) with those of numerical diagonalisations of fi- 
nite matrices; the agreement is good already for N = 32. 

We have also found a way to calculate exactly the dis- 
tribution /(r) of the width of the most narrow resonance 
among those falling in a window [E — W/2, E + W/2] in 
the vicinity of a given point E in the spectrum. Assum- 
ing that the mean number n — W/A of resonances is 
large (n >• 1), but still ^ 1 to preserve spectral local- 
ity (the density of states should not change significantly 
across the spectral window): 



/(r) 



A 



(8) 



This distribution is of great interest in the theory of ran- 
dom lasingB. The functional form of the distribution was 
found in Ref. || by employing plausible qualitative argu- 
ments yielding Eq. (j^), but with renormalised effective 
coupling g replaced by its "weak coupling" limit ^ jlixv. 
We see that the difference with exact formula amounts 
to the factor 2 in the exponent for the case of perfect 
couphng 7 = 1. In Fig. ||, the result (||) is compared 
to results of numerical diagonalisations for = 128 and 
W = 0.2A, corresponding to n « 8.15. 

In the remainder of this article, we outline the deriva- 
tion of the results (|),(0),(|). The main idea is to use 
that fact that the complex eigenvalues (resonances) 5^ 



are poles of the M x M scattering matrix S(£') in the com- 
plex energy plane. Using the standard expression for the 
scattering matrix in terms of the non-Hermitian Hamil- 
tonian Hat (see e.g. Ref. ^) the residues corresponding 
to these poles can be found and after some algebraic ma- 
nipulations we arrive at the following relation: 



Tr <^ Res 



S{£) 



Res 



(9) 



(^m ^n) (£n ^m) 



This relation is valid for arbitrary M, but for M > 1 
it appears to be of no obvious utility, due to difficulties 
in evaluating the ensemble average of the trace of the 
residues on the left-hand side. However for the case of one 
single open channel the scattering matrix can be written 
as 



N 



si£) = n 



k=l 



£~£k 



N 



k=l 



£* - £t 



(10) 



which follows, up to an irrelevant "non-resonant" phase 
factor, from the requirement of S~ matrix analyticity 
in the upper half-plane and unitarity for real energies. 
Substituting Eq. ( [lo| ) into Eq. (jp) yields the relation: 



{£n £„) i^m £rn) JJ^ £k JJ^ 



-£k 



{£n £*n) 



f _P, 11 f*^f* 
k^n ^" k^m ^™ 



(11) 



expressing the eigenvector overlap matrix in terms of 
complex eigenvalues £k^- This gives a possibility to find 
the diagonal and off-diagonal correlators, Eqs. (||,^), by 
averaging Omn overJinown joint probability density of 
complex eigenvaluesE3 for the single-channel scattering 
system: 

r{£u-,£N) = ^p^|A{£i,...,£^}|2 (12) 
fc 

Using this expression one may notice that 

0{£) = 21^e-M^^7^-(^^-i)7i]e-f (13) 



x(det(f-7^t)(5*_^))^^_^ 

where TYat-i stands for the non-Hermitian matrix 7i of 
the same type a.s Hn but of the lesser size {N— 1) x {N — 
1), and with coupling 7 replaced by a modified coupling 
71 = 7 — Im £ . Analogously 



0(^1,^2) 



~ JV-3 

72 „_ir7V72-(Ar-2)72l 
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(14) 



xe-f5:^^.(^S+^;^)(£i-£r)(f2-f2*) 

X (det (£1 - n^) (£* - n^) {£2 - n) {£; - n))^^ 
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where Hn-2 is of the size (iV — 2) x {N — 2), and with 
coupHng 7 replaced by a modified coupling 72 = 7 — 
Im £1 — Im £2 ■ 

In this way the problem is reduced to calculating a 
correlation function of characteristic polynomials of large 
non-Hcrmitian matrices. A closely related object was 
calculated in Ref. |lO[ and we can adopt those methods 
to our case. The scaling limit ^ 1 such that Imfi_2 ~ 
Vx.2 ^ 2VL = Re(£i -£2) A oc N^^ of the resulting 
expressions yields the formulas Eqs. (||)-(0) above. 

Let us briefly comment on a way of calculating the 
distribution Eq. (^) of the widths of the most narrow 
resonance in a given window. Instead of extracting such 
a quantity from the joint probability density Eq. ( |l2| ) we 
find it more convenient to consider 

7'(zi,...,z„)cx ^|A{zi,...,z„}|2 (15) 

n 

x<5(l-T - []lm|zfe| 
fc=i 



defined for complex variables = T'ie^* inside the unit 
circle: — \zi\ < 1. For < T < 1 this formula has 
interpretation of the joint probability density of complex 
eigenvalues Zi for the ensemble of n x n subunitary matri- 
ces and is a very natural "circular" analogue of Eq. (12). 
The similarity is in no way a superficiaLane, but rather 
has deep roots in the theory of scatteringtZl. The parame- 
ter T controls the deviation of the corresponding matrices 
from unitarity, much in the same way as the parameter 
7 controls the deviation oi TLn from Hermiticity. More 
precisely, T should be associated with the renormalised 
coupling constant g via the relation g = 2/T — 1. In the 
limit n 3> 1 the eigenvalues Zk are situated in a narrow 



vicinity of the unit circle. Their statistics is shown0 to 
be indistinguishable from that of the complex eigenen- 
ergies 5, when the latter considered locally, i.e. on the 
distances comparable with the mean spacing A. In par- 
ticular, the distances 1 — from the unit circle should 
be interpreted as the widths of the resonances. 

The form of Eq. ( [l5| ) allows one to integrate out the 
phases Oi by noticing that: 



27r 



rue 



k<j 



{a} 



„2a„ 



where the summation goes over all possible permutations 
{a} — {ai, an) of the set l,...,rt (in fact in the right- 
hand side we deal with the object known as "permanent" , 
see e.g. Ref. ^|). In this way we arrive at a joint proba- 
bility density of the radial coordinates only, and the dis- 
tribution Eq. (|^) follows after a number of integrations 
and the limiting procedure n S> 1. 

In conclusion, we presented a detailed analytical and 
numerical investigation of statistics of resonances and as- 
sociated bi- orthogonal eigenfunctions in a random ma- 
trix model of single channel chaotic scattering with bro- 
ken time-reversal invariance. Among challenging prob- 
lems deserving future research we would like to mention 
extending our results to the case of more than one chan- 
nel and to time-reversal invariant systemsEEl, as well as 
the problem of understanding fluctuations of the non- 
orthogonality overlap matrix Omn- 
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